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Представлены результаты теоретического исследования процесса теплопереноса при выпадении на поверхность элементов 
конструкции газового тракта котлов частиц летучей золы и результаты численного анализа уровня термических напряжений, 
возникающих в тонком приповерхностном слое металла в малой окрестности частицы. Задача решена методом конечных раз- 
ностей. Установлено, что значения упругих напряжений, возникающих в тонком, прилегающем к поверхности нагрева слое ме- 
талла, достаточно велики. Возможна реализация механизма локальной эрозии металла элементов газового тракта котлов за счет 
растрескивания поверхности с образованием сетки мелких трещин, как в радиальном так и в окружном направлениях, даже при 
инерционном осаждении частиц золы с малыми скоростями движения. 


Известно, что работа современных котельных аг- 
регатов сопровождается сложными, взаимосвязан- 
ными физическими и химическими процессами в 
тракте дымовых газов [1]. Эти процессы оказывают 
сильное влияние на работоспособность металла, из 
которого изготовлены узлы и детали котельного обо- 
рудования [1-3]. Наиболее часто неполадки и пов- 
реждения происходят на водяных экономайзерах, 
топочных экранах, пароперегревателях [2]. Износ 
поверхности труб экономайзеров принято считать 
[2] следствием ударов частиц летучей золы о поверх- 
ность труб. Для снижения масштабов эрозии метал- 
ла труб вводятся ограничения на скорости движения 
дымовых газов. Но при малых скоростях дымовых 
газов возникают золовые заносы, вызывающие рост 
сопротивления и ухудшение теплообмена [1,2]. 

В то же время детальный механизм эрозии ме- 
таллов в рассматриваемых условиях котельного 
оборудования до настоящего времени не установ- 
лен однозначно. В данной работе проведен теоре- 
тический анализ условий реализации одного из 
возможных механизмов эрозии металла элементов 
конструкции тракта дымовых газов, например, 
трубной доски воздухонагревателя. 

Рассматривается механизм эрозии металла, явля- 
ющийся следствием не ударов частиц летучей золы, 
движущихся с большой скоростью под определен- 
ными углами атаки по отношению к поверхности 
доски, а следствием инерционного осаждения таких 


частиц с малой скоростью на поверхность элементов 
газового тракта. Частицы золы всегда имеют сущест- 
венно более высокую температуру, чем, теплообмен- 
ная поверхность, которая охлаждается за счёт отвода 
тепла из приповерхностных "горячих" слоёв металла 
в глубинные слои, имеющие всегда относительно 
низкую температуру. При попадании частиц, нагре- 
тых до температур на несколько сотен градусов вы- 
ше, чем поверхность, на которую частица осаждает- 
ся, в тонком приповерхностном слое металла возни- 
кают термические напряжения, обусловленные гра- 
диентами температуры вблизи поверхности нагрева. 
Эти напряжения возникают в малой окрестности 
частицы и не приводили бы к значительным дефор- 
мациям металла, если бы случаи осаждения частиц 
были единичны. В реальных же условиях при любых 
скоростях движения дымовых газов выпадения час- 
тиц происходят непрерывно. Затем зола может уно- 
ситься потоком или оставаться на поверхности эле- 
мента газового тракта. При этом одиночный случай 
осаждения частицы в определенных условиях может 
приводить лишь к возникновению одиночных ло- 
кальных микротрещин, но частые последовательные 
повторения влекут за собой разрушение металла в 
результате его малоцикловой усталости [4]. 

Цель данной работы - теоретический анализ 
полей температур и упругих напряжений в окрест- 
ности нагретой до высоких температур частицы ле- 
тучей золы, выпавшей на поверхность металла. 
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Постановка задачи 

При решении задачи предполагалось, что час- 
тица имеет форму цилиндра малой высоты. Задача 
решалась в цилиндрической системе координат, 
начало которой совпадало с точкой пересечения 
оси симметрии и внешней границы частицы. Была 
принята следующая схема процесса. Частица выпа- 
дает на поверхность пластины металла с начальной 
температурой Т„ из дымовых газов с температурой 
Т г (при этом Т 0 <Т г ). В пластине выделяется ци- 
линдр, характерная высота которого Н и радиус г г 
много больше соответственно высоты и радиально- 
го размера частицы. Параметры внешнего газового 
потока (коэффициент теплоотдачи и Т г ) известны и 
постоянны на всех границах частицы и металла в 
любой момент времени. Предполагалось, что теп- 
лообмен между высокотемпературной средой и 
частицей (а также между средой и металлом) осу- 
ществляется за счет конвекции и излучения. 

Приняты следующие основные допущения, не 
накладывающие серьезных ограничений на общ- 
ность постановки задачи. 

1. Форма частицы после осаждения не изменяет- 
ся. Для частиц золы, находящихся в момент 
осаждения в твердом состоянии и имеющих от- 
носительно невысокие скорости движения, та- 
кое допущение является обоснованным. 

2. Теплообмен между частицей и металлом осуще- 
ствляется в условиях идеального контакта. 

3. Прилегающий к поверхности слой металла на- 
ходится в твердом состоянии, и после выпаде- 
ния частицы не происходит деформации этой 
поверхности. 

4. Условия теплообмена на всей поверхности час- 
тицы и поверхности металла, остающейся отк- 
рытой, принимаются адекватными. 

5. Теплофизические свойства материала частицы 
и металла не зависят от температуры. 

6. Не учитываются возможные химические реак- 
ции между частицей и металлом. 

Математическая модель 

Система нестационарных дифференциальных 
уравнений теплопроводности для системы "частица 
- металл" (рис. 1), соответствующая сформулиро- 
ванной физической модели имеет следующий вид: 



Рис. 1 . Схема области решения задачи 
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При постановке задачи принято, что оси сим- 
метрии частицы и выделяемого в пластине металла 
цилиндра совпадают. Граничные и начальные усло- 
вия для сформулированной задачи имеют вид: 
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О < 2 <к, г = г х 
дТ 

-А, а ; =а(Т,-Т х ) + 08 пр (Т 4 - Т х 4 ), (10) 

1 = 0, Т х (г, 2 ) = Т 10 , Т 2 (г,г) = Т 0 . (11) 

Значения тангенциальных а, и радиальных ст г 
термических напряжений в приповерхностном 
слое металла рассчитывались по приближенным 
выражениям [5], полученным для условий воздей- 
ствия на поверхность пластины теплового источ- 
ника заданной интенсивности радиусом г х . 

а,.=-[(2ЕР(Т-Т 0 )(г х /г) 3 ]/3(1- И ), г>г х , (12) 


а т =[(ЕР(Т-Т 0 )(г х /г) 3 ]/3(\- И ), г>г х , (13) 


а г =< 7 1 =-[(2Е/3(Т-Т 0 )]/3(\- /л), г<г х . (14) 

Здесь Е - модуль упругости: /3 - температурный 
коэффициент линейного расширения: /и - коэф- 
фициент Пуассона. 


107 




Известия Томского политехнического университета. 2004. Т. 307. № 5 


Метод решения 

Сформулированная система дифференциаль- 
ных уравнений с соответствующими краевыми и 
начальными условиями (1-14) решена методом ко- 
нечных разностей [6]. Для решения разностных 
аналогов двумерных уравнений использовался ло- 
кально-одномерный метод в сочетании с методом 
итераций. Одномерные разностные уравнения ре- 
шены методом прогонки с использованием неяв- 
ной четырехточечной разностной схемы [6]. Сеточ- 
ные параметры выбирались таким образом, чтобы 
для частиц любого размера обеспечить не менее 
50 узлов разностной сетки по каждой координате. 

Численные исследования проведены при следу- 
ющих значениях параметров и постоянных [7, 8]: 
7)=1000 К; Г 0 =400, 500, 600, 700 К; а=200, 500, 
1000 Вт/м 2 -К; й=50, 100, 150, 200 мкм; Г 10 =1000 К; 
Я 2 =46 Вт/м-К; с 2 =500 Дж/кг-К; р 2 =7800 кг/м 3 ; 
А|=1,5 Вт/м-К; ^=1100 Дж/кг-К; р 2 =1700 м/м 3 ; 
І7=2- 10 11 Па: /3= 15- ІО- 6 1/К; р= 0,5. 

Основные расчеты проведены для частиц, высота 
которых равна диаметру. Такое соотношение продоль- 
ного и поперечного размеров характерно для частиц, 
образующихся при сгорании типичных углей. 

Результаты численных исследований и обсуждение 

На рис. 2 приведены типичные результаты рас- 
четов в виде зависимости температуры зоны кон- 
такта Т К = Т(0,Н,1) на оси симметрии от времени для 
частиц золы разного диаметра. На этом же рисунке 
приведена зависимость температуры поверхности 
металла Т п от времени в условиях отсутствия частиц 
при адекватных условиях теплообмена. Видно, что 
выпадение частиц в условиях идеального теплового 
контакта приводит к скачкообразному росту Т к от 
начального значения Г 0 = 400 К до величины «700 К 
(при Т г = 1000 К). Затем, с ростом 1 значение Т к 
уменьшается. Также хорошо видно, что уменьше- 
ние размеров частицы приводит к существенному 
снижению величины Т к в каждый момент времени. 
Этот результат обусловлен, скорее всего, различ- 
ным теплосодержанием частиц в момент выпаде- 
ния. Мелкие частицы "остывают" быстрее крупных, 
отдавая тепло в металл, и одновременно закрывая 
участок поверхности, (приближенно равный по 
площади поперечному сечению частицы), от воз- 
действия внешней высокотемпературной среды. 

Зависимости температуры поверхности частицы 
золы в точке с координатами г = 0иг=0от времени 
7] = Г(0,0,/) приведены на рис. 3 также для их разных 
характерных размеров. Видно, что перепад темпера- 
тур по частице на расстоянии Н значителен для час- 
тиц любого размера (от 100 до 700 К). При этом 
Г/О, 0,/) при И = 50 мкм за 2 мс падает от начального 
значения 700 до 400 К. Температура же крупных час- 
тиц Г/0, 0,7) за этот же промежуток времени даже 
растет почти на 100 К. Данный результат дает осно- 
вание сделать вывод о том, что мелкие частицы игра- 
ют в большей степени роль теплоизолятора, чем ис- 
точника нагрева. Большие же нагревают металл. 



Рис. 2. Зависимости температуры зоны контакта Т к =Т(0,Іт,1) 
на оси симметрии от времени для частиц золы разно- 
го диаметра г, при Т 0 =500 К: 1) 200; 2) 150 ; 3) 100; 
4) 50 мкм; 5) граничные условия 1-го рода 



Рис. 3. Зависимости температуры поверхности частицы золы 
Т,=Т(0,0,і) от времени при Т 0 =5ОО К для их разных 
характерных размеров г,: 1) 50; 2) 100; 3) 150; 
4) 200 мкм 

Установлено также, что пятикратное увеличе- 
ние а приводит к очень незначительным (не более 
3 %) изменениям значений Т к и Г/0,0,/) во всем ди- 
апазоне изменения времени. К аналогичным ре- 
зультатам приводит и увеличение температуры 
внешней среды. Эти результаты позволяют сделать 
вывод, что масштабы влияния частиц на темпера- 
туру зоны контакта определяются в начальный 
промежуток времени только параметрами самой 
частицы (теплофизические характеристики, раз- 
мер) и почти не зависят от параметров внешнего 
высокотемпературного газового потока. 


с г , МПа 



Рис. 4. Радиальные термические напряжения а„ для стали 
при г, =200 мк, где 1) 700; 2) 600; 3) 500; 4) 400 К 


108 




Технические науки 


На рис. 4 приведены типичные распределения по г 
радиальных напряжений о> в сечении і = 0. Видно, что 
численные значения о; достигают 7- ІО 7 Па. При этом 
величина напряжений растет с ростом перепада тем- 
ператур АТ=Т Ч -Т 0 . Но даже при относительно малом 
перепаде АТ( Г,=700 К) величина судостигает 310 7 Па. 
Проводя анализ полученных результатов, можно отме- 
тить, что полученные значения ст г достаточно далеки 
при температурах 300... 500 К от пределов прочности 
<7„ р типичных углеродистых сталей. Но если учесть за- 
висимость ст„ р [9] от температуры (таблица), то можно 
сделать обоснованный вывод о том, что расчетные 
значения а„ р (рис. 4) и табличные величины а пр (табли- 
ца) в диапазоне высоких температур достаточно близ- 
ки. Соответственно, вероятность возникновения в 
тонком приповерхностном слое металла радиальных 
сжимающих напряжений, превышающих предел 
прочности металла при соответствующей температуре 
достаточно велика. При выполнении условия а=а пр 
возможно возникновение в нагретом слое трещин, 
ориентированных по радиальной координате. 


МПа 



Рис. 5. Тангенциальные термические напряжения <т„ для ста- 
ли, при т,= 200 мкм, где 1) 700 ; 2) 600; 3) 500; 4) 400 К 

Таблица. Пределы прочности типичных сталей при высоких 
температурах, о пр , МПа 
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Аналогичные распределения получены и для 
тангенциальных напряжений в прогретом слое 
(рис. 5). Можно сделать вывод о том, что в тонком 
приповерхностном слое стали, например трубной 
доски воздухонагревателя, возможно возникнове- 
ние и окружных термических напряжений, сопос- 
тавимых по величине с пределом прочности стали 
при соответствующей температуре. Окружные нап- 
ряжения могут приводить к возникновению сетки 
трещин (или одиночных трещин) ориентирован- 
ных по окружной координате в малой области вы- 
падения частиц золы. Вероятность реализации та- 
кого механизма эрозии металлов элементов ко- 
тельного оборудования достаточно велика, потому 
что для типичных режимов работы разность темпе- 
ратур дымовых газов (и соответственно частиц ле- 
тучей золы) и поверхности элементов газового 
тракта составляет несколько сот градусов. 


Необходимо отметить, что при проведении числен- 
ного анализа были рассмотрены частицы золы разме- 
рами до 500 мкм. В настоящее время нет оснований го- 
ворить как о доминировании в продуктах сгорания уг- 
ля частиц мелких фракций (с характерными размерами 
единицы мкм), так и о преобладании крупных частиц 
(с размерами до 200 мкм) [10]. Получение достоверных 
экспериментальных данных о массовой концентрации 
и распределении частиц по размерам в продуктах сго- 
рания угля может создать предпосылки для более об- 
щих выводов о масштабах влияния частиц конденси- 
рованной фазы на закономерности процессов прогре- 
ва и эрозии металлов котельного оборудования. 

Заключение 

Численно исследовано локальное температурное 
поле в малой окрестности зоны контакта твердой, 
разогретой до высоких температур частицы, выпав- 
шей на поверхность элемента конструкции котла из 
потока дымовых газов. Полученные результаты ха- 
рактеризуют не только возможный диапазон изме- 
нения температуры зоны контакта для различных 
условий теплообмена с внешней средой, материалов 
и размеров частиц. На основании полученной ин- 
формации можно с достаточно высокой степенью 
надежности прогнозировать реальные режимы 
прогрева металла за счет передачи тепла нагретыми 
до высоких температур твердыми или жидкими час- 
тицами, которые при попадании на поверхность ме- 
талла создают местные центры высоких температур. 

Полученные результаты позволяют также сделать 
заключение о возможности реализации механизма 
локальной эрозии металла в газовых трактах котель- 
ного оборудования. Эрозия может быть обусловлена 
возникновением в тонком приповерхностном слое 
металла термических напряжений (сопоставимых по 
величине с пределами прочности данного металла 
при соответствующих температурах) за счет локаль- 
ного интенсивного прогрева металла выпавшей на 
поверхность элемента конструкции частицы летучей 
золы. Образование в результате реализации такого 
механизма сетки микротрещин на поверхности ме- 
талла влечёт за собой ослабление прочности металла 
в тонком приповерхностном слое, прилегающем к 
поверхности нагрева, и последующее разрушение с 
уносом мелких частиц потоком дымовых газов. 

Список обозначений 

Т - температура, К; с - удельная теплоемкость, Дж/кг-К; 
р - плотность, кг/м 5 ; А - коэффициент теплопроводности, 
Вт/м К; 1 - время, с; г, і - цилиндрические координаты, м; 
/; - высота частицы, м; Н - высота цилиндра выделяемого в 
пластине металла, м; а - коэффициент теплообмена, 
Вт/м 2 ' К; а - постоянная Стефана-Больцмана; Т г - темпера- 
тура внешнего газового потока, К; е„ р - приведенная степень 
черноты; Т 0 - температура металла в начальный момент вре- 
мени, К; Т к - температура в зоне контакта на оси симмет- 
рии, К; 71(0,0,0 - температура поверхности частицы на оси 
симметрии, К; Т п - температура поверхности конденсиро- 
ванного вещества без частицы, К. Индексы 1 и 2 соответ- 
ствуют частице и металлу. 
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